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Multi-orbital Hubbard models are shown to exhibit a spatially isotropic spin-triplet superconduct¬ 
ing phase, where equal-spin electrons in different local orbitals are paired. This superconducting 
state is stabilized in the spin-freezing crossover regime, where local moments emerge in the metal 
phase, and the pairing is substantially assisted by spin anisotropy. The phase diagram features 
a superconducting dome below a non-Fermi liquid metallic region and next to a magnetically or¬ 
dered phase. We suggest that this type of ffuctuating-moment induced superconductivity, which 
is not originating from fluctuations near a quantum critical point, may be realized in spin-triplet 
superconductors such as strontium ruthenates and uranium compounds. 


Spin-triplet superconductivity, in the sense of equal- 
spin pairing, is believed to occur in a number of cor¬ 
related materials. The best candidate is the layered 
compound Sr2Ru04, where the Knight shift remains un¬ 
changed across the superconducting phase boundary, in 
stark contrast with the behavior expected for spin-singlet 
pairing [ij. In the iron pnictides, where a spin-triplet su¬ 
perconducting phase has been proposed in early theoret¬ 
ical works 0, the experimental evidence points toward 
spin-singlet pairing, although in LiFeAs a spin-triplet sce¬ 
nario is still being debated [1, Q. The uranium based 
superconductors are also possible candidates for spin- 
triplet pairing. In compounds such as UGe2, UCoGe 
and URhGe, the superconducting state is found near a 
ferromagnetic phase and the two orders may even coex¬ 
ist @-@1. For a deeper understanding of unconventional 
superconductivity in strongly correlated electron systems 
with multiple active orbitals, it is thus important to clar¬ 
ify the mechanisms which can lead to spin-triplet pairing. 


While a p-wave symmetry is usually assumed for the 
pairing state in spin-triplet superconductors, an s-wave 
spin-triplet pairing is also possible by taking into ac¬ 
count the orbital degrees of freedom. The mechanism of 
this unconventional superconductivity can be easily un¬ 
derstood iiEl: same-spin electrons tend to occupy 
the same site due to the Hund coupling which favors 
high-spin states. A new insight in this paper is that 
the s-wave spin-triplet pairing is closely connected to 
the emergence of local magnetic moments in so-called 
Hund metal s fl^. T his class of materials, which includes 
ruthenates [Ulil and iron pnictides [lol - l^ . exhibits 
Hund-coupling induced correlation effects and character¬ 
istic non-Fermi liquid properties. The underlying phe¬ 
nomenon is spin-freezing [l7| : In a narrowly defined range 
of fillings and interaction strengths, long-lived magnetic 
moments appear in the metal phase of multi-band sys¬ 
tems with Hund coupling (formation of a large composite 
spin, see right panel of Fig.[T](a)). In the absence of long- 
range order, the emerging local moments will be screened 
at sufficiently low temperature, so that there is no quan¬ 
tum phase transition associated with spin freezing. How¬ 


ever, screening large local moments is difficult, and the 
Fermi liquid coherence temperature becomes very low 
[3, li^- Hence, as demonstrated here, a spontaneous 
symmetry breaking pre-empts the screening of the mo¬ 
ments. While deep in the spin-frozen regime the long- 
lived local moments order magnetically at low tempera¬ 
tures, the emerging and fluctuating local moments in the 
spin-freezing crossover regime generate spin-triplet pair¬ 
ing. This leads to the formation of a superconducting 
dome separating the Fermi liquid metal from the mag¬ 
netically ordered region and results in phase diagrams 
which closely resemble those of unconventional supercon¬ 
ductors. 

We consider a three-orbital Hubbard model whose 
Hamiltonian is given by 

= 'y ~ + U ^ ^ ( 1 ) 

Ac 7 £T 27 

+ u' E + {U'-J) Y. (j 

2(T,7<7' 2(7,7<7' 

— aJ ^ + H.C.), 

i,7<7' 

where i is the site index, 7 = 1,2,3 the orbital index, 
a =t, \. the spin index, and a represents the complemen¬ 
tary spin (f =),). Ek is the dispersion of electrons on 
the lattice, and ^ is the chemical potential. The interac¬ 
tion terms contain the intra-orbital iU) and inter-orbital 
{U') Coulomb repulsions, and the Hund coupling J. The 
parameter a controls the anisotropy in spin space, i.e. 
a = 1 corresponds to a spin-rotationally invariant system 
and a = 0 to the Ising anisotropic case where the interac¬ 
tions are only of density-density type. A spin anisotropy 
may originate from spin-orbit coupling, and the param¬ 
eter a allows us to incorporate this effect in a simple 
manner [2^ . 

An anisotropic coupling in spin space should in prin¬ 
ciple only change the prefactor of the spin-flip term in 
Eq. 11]). However, for J > 0 the pair-hopping term, which 
transfers two electrons in the same orbital to another or¬ 
bital, is not important. This is because the pair-hopping 
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FIG. 1: (a) Illustration of possible local configurations with two electrons in three orbitals, and the corresponding energies 

E. (b) Filling dependence of the local magnetic susceptibilities for U = 1 and T = 0.005. The black curve shows the local 
magnetic susceptibility Xioc and the blue curve the contribution from long-lived (frozen) moments, (c) Filling dependence of 
the renormalization factor 2 for the same parameters. 


favors the state shown in the left panel of Fig.[lja), which 
is hardly realized due to the presence of the intra-site 
Coulomb interaction U which is larger than J. We there¬ 
fore consider it more convenient to put the anisotropy 
factor in front of both terms, so that a interpolates be¬ 
tween the familiar Ising and rotationally invariant limits. 

While the phenomena discussed in this paper are 
generic features of multi-orbital systems with non-zero 
Hund coupling parameter J, we will show results for the 
3-orbital case with U' = U — 2J, J/U = 1/4 and consider 
a semi-circular density of states with bandwidth W = 1. 
(We neglect specific material effects related to the partic¬ 
ular shape of the density of states and to the hybridiza¬ 
tion between different orbitals, aiming at a general, ma¬ 
terial unspecific, discussion of the physics.) The model 
is solved using the dynamical mean field theory (DMFT) 
, combined with a numerically exact continuous-time 
Monte Carlo method [1^ [1^ . This formalism captures 
local correlation effects. 

To illustrate the spin-freezing phenomenon, we com¬ 
pute the dynamic contribution to the local magnetic sus¬ 
ceptibility 

rP 

Axioc = / dr ((5 ,(t)5,) - (5,(/3/2),5,)), (2) 

Jo 

where /3 = 1/T is the inverse temperature and ^(r) = 
e^^^e”^-^. The operator Si = ^ ~ 

M = 3 the number of orbitals, mea¬ 
sures the local spin. The first term on the right-hand 
side yields the local magnetic susceptibility (xioc)- In 
Eq. (m we subtract the long-time correlator {Si{/3/2)Si), 
which reflects the magnitude of long-lived frozen mo¬ 
ments [i3,|3^. Hence, the quantity Ayioc measures the 
fluctuations of the moments. Figure [ijb) shows the fill¬ 
ing dependence of these quantities for a = 0. While 
the local susceptibility xioc monotonically increases with 
increasing n, the fluctuation Axioc reaches a maximum 
at n ~ 1.9. This peak indicates the crossover between 


the Fermi liquid and spin-frozen regimes, and we use the 
location of the maximum as our definition of the “spin¬ 
freezing crossover.” The spin-freezing is also reflected in 
the renormalization factor z, or mass-enhancement fac¬ 
tor l/z, of the quasi-particles [l^. For the estimation, 
we use the ansatz S(a; 0) = a + buj and determine the 
coefficients by fitting the numerical data. Specifically, 
we fit the self energy by the Fade approximation using 
the lowest two Matsubara frequencies, and compute the 
renormalization factor by the relation z = (1 — 6)“^. Fig- 
ure[ljc) exhibits a drop of 2 in the spin-freezing crossover 
region. Qualitatively similar results are obtained in the 
SU(2) symmetric case {a = 1, dashed lines). 

To study the stability regions of ordered phases, we 
calculated the susceptibilities 



where N is the total number of sites. The operator is 
given by 

r (FM) 

^={ (AFM) (4) 

for 7 7' (SC) 

for ferromagnetic order (FM), anti-ferromagnetic order 
(AFM) and s-wave inter-orbital spin-triplet superconduc¬ 
tivity (SC). Xi is a sign which depends on the sublat¬ 
tice. A divergence in Xff (or equivalently a sign-change 
in 1/x^) indicates a possible transition into a long-range 
ordered phase. The susceptibilities in Eq. (|3]) can be de¬ 
rived from the two-particle Green’s function [i^. We 
calculate the vertex part of this Green’s function from 
the local impurity problem, and obtain the lattice two- 
particle Green’s function by solving the Bethe-Salpeter 
equation. While we have also calculated the susceptibil¬ 
ities for other types of orders, such as orbital ordering, 
only the quantities listed in Eq. ((T|) diverge in the param¬ 
eter regions considered in this paper. 
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FIG. 2: (a) Filling dependence of the inverse susceptibilities for FM, AFM and SC orders for U = 0.75, T = 0.005 in the 

system with Ising spin anisotropy (a = 0). A negative 1/xo indicates a long-range ordered phase with order parameter O. 
(b,c) Interaction-filling phase diagrams at T — 0.005 and T = 0.0025, respectively. The black dashed line shows the location of 
the spin-freezing crossover in the system without long-range order. 


Figure [2ja) shows the inverse susceptibilities for T = 
0.005, U = 0.75, a = 0 and different fillings. Symme¬ 
try broken phases exist in the regions where 1/x^ < 0. 
Repeating this analysis for different U, we obtain the 
T = 0.005 phase diagram shown in Fig. [2Ib). At 
U > 1.25, a FM phase appears 3l|, while near half-filling 
the AFM phase is stable. A new result is the existence of 
a SC region connecting the FM and AFM phases. This 
spin-triplet SC phase is clearly associated with the spin¬ 
freezing crossover, indicated by the black dashed line, 
which suggests that the fluctuating local moments at the 
border of the spin-frozen regime induce the pairing. We 
also show the phase diagram at a lower T = 0.0025 in 
Fig. He), where the SC region expands. 

We note that there is no direct attraction among elec¬ 
trons, although the superconductivity is realized by form¬ 
ing electron pairs. The effective attraction for Cooper 
pairs can be understood from the imbalance of the 
Coulomb interactions 3 32-3^. To see this, let us con¬ 
sider the situation with two electrons on one lattice site. 
As shown in Fig. there are three kinds of config¬ 

urations. The energetically most favorable one is the 
electron pair with the same spins, because of the Hund 
coupling J. Hence, two same spin electrons will tend to 
occupy the same site (but different orbitals), even though 
a repulsive Coulomb interaction acts between them. In¬ 
deed, the local effective interaction U among same-spin 
electrons can be derived within second-order perturba¬ 
tion theory [3], which results in 


U^{U' -J)- [2UU' + {V - Jf + U'^\x\oc, (5) 


with xioc = Ayioc in the weak-coupling approximation 
[diij l. If the second term on the right-hand side dominates 
the first-order term due to a large xioc the interaction 
becomes attractive. 


A reduction of the interaction energy below 17 — 3 J is 
not possible if the number of electrons is constrained to 
two per site, on average. On the other hand, the electrons 
do not have to occupy the same site when n < 1, and the 
superconducting state is never realized in this case. 

Next, we discuss the temperature-filling phase diagram 
shown for U = 0.75 in Fig. Ha). With hole doping from 
half-filling (n = 3), the AFM transition temperature de¬ 
creases and becomes zero at n ~ 2.3, which is close to 
the spin-freezing line. By further doping with holes, we 
find the spin-triplet SC phase with a dome-shaped T^. 
If we fix the filling to n = 2 and change C/, we obtain 
the phase diagram shown in Fig. Hb). With decreasing 
U (increasing pressure, experimentally), the FM order is 
destroyed and again a SC dome appears next to the mag¬ 
netic region. The light blue diamonds in Fig. [3] show the 
points where the inverse pairing susceptibility reaches its 
minimum as indicated by the arrow in Fig. H^-)- Above 
Tc, this corresponds to the maximal pairing instability, 
while below Tc this minimum approximately locates the 
maximum order parameter (4fl| . The close resemblance 
with the spin-freezing crossover line shows the relevance 
of spin freezing for the present superconductivity. 

SC domes are usually understood as a manifestation 
of fluctuations associated with magnetic quantum crit¬ 
ical points. However, the superconductivity revealed in 
this paper is induced by local magnetic fluctuations in the 
spin-freezing crossover regime. Nevertheless, the SC or¬ 
der naturally appears in the vicinity of a magnetic phase, 
since the strengthening of the magnetic moments deeper 
inside the spin-frozen regime causes magnetic ordering. 
Furthermore, the normal state above the SC dome is a 
non-Fermi liquid whose properties are influenced by the 
spin-freezing crossover |l7l |. 

So far we have shown results for the system with Ising 
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FIG. 3: (a) Temperature-filling phase diagram at t/ = 0.75. 

(b) Temperature-interaction phase diagram at n = 2. The 
black dashed line indicates the spin-freezing crossover in the 
system with suppressed long-range order. The light blue dia¬ 
monds show the minimum of the inverse pairing susceptibility. 
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FIG. 4: Dependence of the superconducting transition tem¬ 
perature on the spin anisotropy a for U — 0.875 and n = 2. 


anisotropy. We now clarify how the superconductivity is 
affected by the spin-flip term in the model with a 0. 
Figure m shows Tc for U = 0.875 and filling n = 2. As 
a is increased from 0, the transition temperature de¬ 
creases and drops below the lowest accessible tempera¬ 
ture at a = 1. The destabilization of the electron pairs 
by the spin-flip term can be intuitively understood by 
looking at Fig. (Ha). Since the spin-flip term exchanges 
t and 4 , spin electrons residing in different orbitals, the 
configuration shown in the middle panel is favored. As 
a result, the probability for the equal-spin state (right 
panel) decreases compared to the spin-anisotropic case. 
Near a = 1, the fluctuations among the 3 degenerate 
spin-triplet states further destabilize the equal-spin pair¬ 
ing. In fact, a previous study of a single-orbital model 
based on the Eliashberg theory demonstrated the im¬ 
portance of longitudinal fluctuations for pairing, which is 
consistent with our results. 

Finally, let us comment on the potential implications 


of these findings for unconventional multi-band super¬ 
conductors. Because an Ising-type spin-freezing is un¬ 
derlying the fluctuating-moment induced spin-triplet su¬ 
perconductivity, it may be expected to occur in electron 
systems with strong spin-orbit coupling. Promising can¬ 
didates are the uranium-based superconductors UGe 2 0 ], 
URhGe 0 and UCoGe Q, which exhibit a SC phase bor¬ 
dering a FM phase. In these compounds a strong Ising 
spin-anisotropy is observed: the magnetization along the 
easy axis is several times larger than along the hard 
axis 0. Futhermore, Ising-type spin-fluctuations are 
important for the superconductivity, since a magnetic 
field along the magnetic moment destabilizes the pairing, 
while it is much more robust against fields perpendicular 
to the moment @,[ 3 ^. The 5/ electrons in the U ions, 
which play a central role in the low-temperature behav¬ 
ior, have a relatively itinerant nature and are strongly 
correlated. Although a realistic description should in¬ 
volve an Anderson lattice, the local interaction is of the 
Slater-Kanamori type, and hence we expect the same 
spin-triplet SC state. We thus believe that our mech¬ 
anism could be realized in these uranium-based super¬ 
conductors. 

Sr 2 Ru 04 is another candidate compound which might 
exhibit a fluctuating-moment induced superconductivity. 
Here, the spin-orbit coupling is nearly 100 meV [l[, and 
as shown in Fig. [4] the spin anisotropy need not be very 
large to realize a SC state at low temperatures. Also, 
the estimated 17 ~ 0.8 [l^ and the filling n = 4 (same as 
n = 2 due to particle-hole symmetry) place this material 
in the parameter regime where the SC state is found near 
a FM phase (Figures [2Kb) and (c)). A related conipound, 
SrRuOa, with a larger U, becomes a ferromagnet [1^ and 
exhibits the non-Fermi liquid behavior associated with 
spin freezing in the high-temperature phase (ttI . . 

In the iron pnictides, the Coulomb interactions and fill¬ 
ings on the Fe sites can also be close to the spin-freezing 
crossover values [l^, [2l[, and for LiFeAs, in particular, 
the experimental signatures fully support this interpreta¬ 
tion [39|. On the other hand, the 3d electrons have weak 
spin-orbit coupling, and hence a small spin anisotropy. 
Thus, for iron pnictides (including LiFeAs), one can ex¬ 
pect other types of pairing driven by e.g. Fermi-surface 
nesting mechanisms to dominate. 

For an experimental detection of the present super¬ 
conductivity, it is necessary to measure both the spin 
and spatial parts of the pairing state. The spin-part can 
be determined by measuring the magnetic susceptibility, 
e.g. with NMR. If it is identified as spin-triplet, one still 
has to distinguish between p-wave pairing and the pro¬ 
posed s-wave inter-orbital pairing. The difference lies in 
the presence/absence of a node in the gap function which 
will be reflected in a power-law/exponential temperature 
dependence of physical quantities. Our work provides a 
general guiding principle in the search for new uncon¬ 
ventional multi-band superconductors, namely the com- 
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bination of emerging local moments in the spin freezing 
crossover regime and spin anisotropy in heavy elements. 
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Effective attraction from purely repulsive interactions 

Here, we explain how local fluctuations can induce a pairing among repulsively interacting electrons. We follow 
Ref. 1 that deals with a three-component fermion system. In the weak-coupling regime, the effective interactions that 
incorporate bubble diagrams can be written in general form as 

Uafiisi) — UoiP ^ ^^ UaaiXai{Q)Uail3{Q)t ( 1 ) 

ai 

where q = (q, Wm) with q denoting the wave vector and Vm = a bosonic Matsubara frequency. For the present 

three orbital Hubbard model, the indices are given by a = ( 7 , a) where 7 = 1, 2, 3 and a ='|', f. The bare interactions 
are given by U^a-jcT = 0, = U, = U', = U' — J (7 ^ 70- The dynamical susceptibility is defined 

by Xa(q) = — 9 a{k)ga(k + q) where ga(k) is the single-particle Green function for electrons with flavor a. For the 
case of degenerate orbitals considered in our paper, we do not need the index a in the susceptibility. 

In the DMFT approximation, only the local part of the vertex corrections is taken into account [2]. Hence we 
replace the susceptibility by the local one, xioc(ii^m)- (This replacement is not essential for the pairing: the effective 
attraction can be derived even when we consider the q-dependent susceptibility, as discussed in Ref. 1.) By solving 
Eq. ID), the static interaction among 1 f and 2 f electrons can be explicitly derived as 

fr ^ U'-J+(J^- 2UU' - 2I/'J)xioc + iU' - J){U^ - 2J^ + AU'J)xL (0) 

) [!_([/_ J)yi„,][l -{U + 2J)xioc][l + iU- 2U' + J)xioc][l + {U + W - 2J)xioc] ’ ^ ’ 

where we consider the static component: xioc = Xioc(O). The diagrams up to second-order in the interactions are 
shown in Fig. [D In this approximation the effective interaction is given by 

t/it-2t(0) ^U'-J- [2UU' + ([/' - + U'^]xioc- (3) 

Thus if the second-order terms dominate the bare interaction U' — J, the effective interaction can become 

attractive even though the bare interaction is repulsive. Hence, Eq. |3|) shows that strong local fluctuations induce a 
pairing among electrons. This argument is valid in the case of weak interactions, where no local moments are formed. 
In this regime, the relation Ayioc = Xioc holds (see Eq. (2) in the main text for the definition of Ayioc)- 

In the above argument, the local susceptibility can be identified as the magnetic and charge susceptibilities, which 
have the same value in the weak-coupling limit. With increasing repulsive Coulomb interactions, the magnetic 
susceptibility is enhanced and the charge one is suppressed. Hence we expect that in the regime considered in the 
main text, the local magnetic fluctuations primarily contribute to the pairing among electrons in the multi-orbital 
Hubbard model. Indeed our DMFT-fCTQMC calculations demonstrate a clear connection between superconductivity 
and local spin susceptibility. We note that the present discussion cannot be applied to the local-moment regime with 
Axioc ^ Xioc- In this case the expansion from the strong-coupling limit should work as an effective theory. 

The other effective interactions can be derived in a similar manner. Figure [2] shows the local-fluctuation (sus¬ 
ceptibility) dependence of the effective interactions for several values of J/U. Note that we have used the relation 
U' = U — 2J. Here we have only the three independent components Ui-f-ii, ^ 172 ^ and Ui^- 2 ^. As shown in Fig.[Da,b), 
for vanishing or small J, the repulsive interaction is reduced by the screening effect from the other orbitals but stays 
positive (repulsive). On the other hand, for larger J, see Fig.[2|c,d), f/ 1714 , and f/ 1724 , are not much modified, while 
strongly reduced with increasing fluctuations and finally becomes attractive. This behavior naturally explains 
the superconducting mechanism in the weak-coupling regime. 



FIG. 1: Effective attractive interactions from bubble diagrams up to second order. 
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FIG. 2: Effective local interactions as a function of the local susceptibility. 


Details of the susceptibility calculations 


Here we explain the details how to calculate the susceptibility in the framework of DMFT [2]. For diagonal and 
offdiagonal orders, the relevant two-particle Green functions are given by 

xfr7%V'('ri:T-2,'r3,r4) = {Trcl^^{Ti)c^jaiT2)cly^,{T3)CjYa'{T4)) - ((ti (t 2 )) (Frc]y (t3)CjV<t'(T 4 )) , (4) 

X°^jaYa'{Tl,T2,T3,T4) = (^r (ti )c|y (t2)CjV <t'(t 3 (r4)) , (5) 

respectively. Here T,- represents the imaginary time-ordering operator, and we subtract the disconnected part in 
Eq. (U) . We also define the Fourier transformation with respect to time by 


X?j,-yaYa'(.^^n,i£n') = ^ dTidT2dT3dT4 X?j,-ya-y'a'{^1, ^2, T3, T^) 


( 6 ) 


where a denotes ‘diag’ or ‘offd’. Since we are interested only in static susceptibilities, we set the bosonic frequency 
to zero. 

Now we employ the Bethe-Salpeter equation which relates the two-particle Green function and vertex parts. Since 
the vertex part F is local in DMFT [2], it can be evaluated from the local two-particle Green function as 

Xz 2 , 7 cr 7 ^CT') Xz 2 , 7 ct 7 ^{t' ^ ^ ^ ^ X 22 , 7 cr 7icri )^z,71(7172^2 )Xz2,72(727^(7^ (^^^2 * 


^l7icri 71272^^2 


(7) 


Here the two-particle Green function without vertex parts is written as The site index can be neglected if the 

system is uniform. The local quantities are calculated from the effective impurity system using the standard technique 
of GTQMG [3]. Equation ([7]) may be regarded as a matrix equation with respect to the set (n, 7 , cr), and by solving 
this the vertex parts are obtained. In our calculation, the number of Matsubara frequencies for this two-particle 
quantity is typically taken as 150. Although the value of the susceptibility varies as a function of this cutoff, the 
divergent points (phase boundaries) are insensitive to it. 

We note that additional contributions to the vertex parts, which cannot be written in the form of Eq. o, must be 
considered when we have e.g. offdiagonal hybridizations. In the present model, however, Eq. o is sufficient for the 
derivation of the susceptibilities listed in Eq. (4) in the main text. 




















The local vertex extracted above is inserted into the non-local Bethe-Salpeter equation 

k ni7l(7l 71272(72 


y'a' (i^Ti 


, 1£7 


( 8 ) 


In practice, it is convenient to perform the Fourier transformation with respect to the site index before solving the 
matrix equation. Thus we derive the two-particle Green function, which contain the information of the susceptibilities. 
For example, the ferromagnetic susceptibility is given by 


^ M ^ ^ ^ ^' 

ij nn' 77 ' <7<7' 


aa^ a'a 


fie 

A'ij,7cr7'(j' 




(9) 


where is a z-component of spin-1/2 Pauli matrix. On the other hand, the inter-orbital/spin-triplet pairing suscep¬ 
tibility is given by 

^sc = ^ X?®7t7't(i^"’ is„(), (10) 

iji nn' 


with 7 ^ 7 '. The original definitions are given by Eqs. (3) and (4) in the main text. The other susceptiblities for e.g. 
orbital order can be calculated in a similar manner. 


Interpretation of the minimum in the pairing susceptibility below Tc 


Let us consider the meaning of the negative pairing susceptibility obtained in certain filling or interaction regimes. 
According to Landau theory, the free energy is given by 




( 11 ) 


where is the pairing amplitude and b > 0. Since we have the relation 

“ = Xsc/2 


( 12 ) 


between the coefficient a and the pairing susceptibility xsC) the sign reversal of a signals a second-order phase 
transition as shown in Fig. [3] If a < 0, the pair amplitude is given by 


= 



(13) 


This indicates that the magnitude of the order parameter can be deduced from the value of a by assuming that 6 is a 
constant, even though we do not enter the symmetry broken phase. More specifically, the minimum of the (negative) 
inverse susceptibility coincides with the maximum of the pairing amplitude. Although this argument is valid 
only near the transition temperature, we expect that the connection between minimum Xsc ®'^'l maximum pairing 
amplitude approximately holds in a wider temperature range. 



FIG. 3: Landau free energy landscape above and below the transition temperature. 
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Additional data for the phase diagram 


In the main text, we discuss the superconductivity by choosing several interaction parameters. Here we show that 
the behavior remains qualitatively the same when we change the parameters. Figured shows the filling and anisotropy 
parameter a dependences of the inverse susceptibilities and phase diagrams at C/ = 1. The results are qualitatively 
similar to those for U = 0.75 and U = 0.875 shown in the main text. 
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0 0.2 0.4 0.6 0.8 1 
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FIG. 4: Results for ?7 = 1: (Left) filling dependence of susceptibilities at T = 0.005, (Middle) filling-temperature phase 
diagram, and (Right) spin anisotropy (a) dependence of the superconducting transition temperature at n = 2. 
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